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ABSTRACT 


A linear shif t-invarieint image preprocessing technique is 
examined which requires no specific knowledge of any parameter 
of the original image and which is sufficiently general to allow 
the effective radius of the composite imaging system to be arbi- 
trarily shaped and reduced, subject primarily to the noise power 
constraint. In addition, the size of the point-spread function 
of the preprocessing filter can be arbitrarily controlled, thus 
minimizing truncation errors. 



I . INTRODUCTION 


The general problem of image processing has received much 
attention within the last decade. The intense interest , in this 
area arises from the need for the highest possible image quality 
in the increasing application of the many forms of imagery, from 
X~rays in medicine to data collected from satellite-based multi- 
spectral optical line scanners for monitoring earth. resources, to 
the solution of various related problems in many fields of science 
and engineering , made feasible by recent improvements in digital 
computer hardware. The general area of image processing may be 
divided into three major categories; image preprocessing, effi- 
cient image coding, and pattern recognition. There are several 
comprehensive tutorial surveys which cite the significemt tech- 
niques for handling problems in each of these categories and 
which contain extensive references. (1, 12, 23) 

Since no image collecting or imaging system will produce a 
perfect replica of the original image, some further processing is 
usually required. Image preprocessing deals primarily with the 
problem of processing the output of an imaging system in such a 
way that the significant parameters or features of the original 
image are, in some sense, enhanced or restored. This processing 
may be linear or nonlinear, shift-variant or invariant depending 
upon the type of degradation produced by the imaging system. 
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II. STATEMENT OF PROBLEM 


The purpose of this research is the development of a tech- 
nique to reduce the effective aperture radius of multispectral 
optical line scanners used for remote sensing of earth resources. 
The data gathered from such systems is used principally for the 
classification of individual resolution elements by pattern recog- 
nition techniques. The accuracy of such classification techniques 
is usually based upon the supposition that each resolution element 
-of the imaging system output exactly represents a sample of a 
correspondingly located element of the original image. Because 
of the finite aperture size of the scanner, which is not only a 

function of the optics of the scanner but also the impulse response 

(14) 

of any analog signal conditioning or recording equipment, a 

two-dimensional spatial smearing or blurring of the original image 
is produced. This type of imaging degradation essentially maps 
many points from the original image into a single resolution ele- 
ment. In other words, a single resolution element of the imaging 
system output represents a two-dimensional weighted sum of many 
points adjacent to the correspondingly located sample of the 
original image. Thus# depending upon the density and shape of 
the aperture and the spatial and multispectral characteristics 
of the original image, serious classification errors may result. 

This smearing has been observed to seriously affect classification 
accuracy within several aperature diameters of the boundaries of 
data classes. In addition, the classification accuracy of any 
topographical feature of approximately two aperture diameters or 
less; for example, roads, streaims, and buildings at an altitude 
of 1.5 kilometers or more, is substantially reduced. 

It is expected that a reduction of the aperture radius will 
decouple the spatial correlation between adjacent resolution ele- 
ments of the imaging system output, thus correcting each resolution 
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element, so that it more accurately represents a single sample of 
a correspondingly located element of the original image and not 
a weighted sum of adjacent points. Consequently classification 
accuracy of both small topographical features and in the bound- 
ary areas of large data classes should be improved, as well as 
overall spatial resolution of the imaging system output. 

(14 

An analysis of the multispectral optical line scanner system 

indicates that the imaging system degradation could be assumed to 

(20 23 ) 

be linear shift-invariant. ' The proposed preprocessing 

technique is based upon this assumption. The principle advantage 
in making such an assumption was to reduce the cost of preprocessing. 
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III. COMPARISON OF PROPOSED TECHNIQUE TO EXISTING TECHNIQUES 

Numerous techniques have been proposed for 
10 , 15, 19, 22, shift-invariant degraded 

images. The majority of these techniques require some knowledge 

( 23 ) 

of the original image. For exeunple, when the mean-square 

error of the processed image is minimized, which incidently is 

(23) 

not a very effective performance criterion, the resulting 

filter requires a knowledge of the power spectral density of the 
original image. ’ ' ' ' ' When precise knowledge of 

the required parameter of the input signal is not known, the 
resulting error produced by the processor may often negate any 
possible image improvement. Such techniques must necessarily 
require that a different processing filter be used for each 
specific image class comprising images having similar "a priori ” 
statistics. 

In view of the potentially large number of image classes 
comprising the data processed at LARS, the cost of such a pre- 
processing technique requiring a separate "matched” filter for 
each specific image class would be prohibitive. The technique 
examined in this research does not require specific information 
about the original image. Thus a single processing filter for 
all image classes would be required. However, it should be noted 
that the resulting processing filter is suboptimal in the sense 
that "a priori " statistics of the original image are ignored. 

The fundamental objectives of this technique are similar to 

(24) (25) 

those examined by Smith and Stuller. The diagram of 

the basic image preprocessing system is shown in Figure 1. The 
problem is to determine the optimal preprocessing filter point- 
spread function, h^(v), which will make the composite imaging 
system point-spread function g(v), arbitrarily close to an impulse 
function, subject to a constraint On the meah-square noise compon- 
ent in the processed image, n^(v). 
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Stated more precisely, the problem is to choose the h (v) 
that will minimize the functional 

g^(v)dv where v is a two-dimensional vector 

emd 

gR = hj^(v)*hj.(v) 
subject to the constraints 

00 . 

^2 * I 9^(v)dv. 

.00 

The function w(v) is a penalty function designed to force the 
composite imaging point-spread function , g (v) , to be arbitrarily 
duration limited, thus approximating the desired impulse function. 
The more rapidly w(v) increases with increasing v, the more 
rapidly g(v) will decrease with increasing v. Both Smith and 
Stuller chose w(v) = v^, because of resulting mathematical conven- 
iences; although a more general formulation allowing for a higher 
order penalty function would be desirable. It would provide the 
filter designer with an additional parameter for controlling the 
degree of resolution improvement. 

In practice it is also desirable to have h^(v) duration limi- 
ted. Ultimately ciny preprocessing will be performed digitally; 

and since only a finite record length of h (v) may be used, serious 

(12 21 ) ^ 

truncation errors may result. ' The technique proposed 

by Smith did not provide a means for arbitrarily controlling the 
duration of h^(v). The lack of such a constraint also leads to a 
difficulty in obtaining h^(v) from the solution of a differential 
equation. The technique proposed by Stuller provided for an 
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Figure 1 Block Diagram of Imaging and Preprocessing Systems 
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arbitrary control on the duration of h^(v) by allowing the 
solution for h^(v) to contain only a specified number of data 
points. 

The technique proposed in this research adds ein additional 
constraint to the previous two constraints 


K3 = 


s (v)h* (v)dv 
r 


*00 


where s (v) is an arbitrary penalty function designed to duration 
limit the preprocessing filter point-spread function, h^ (v) . The 
addition of this constraint provides a control on the rate of 
decay as well as the duration of h^(v). In addition the proposed 
technique allows for an arbitrary w(v) and obtains a solution 
for h^ (v) by using a different approach from that of Smith or 
Stuller, which may be easily adapted to additional constraints. 
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IV. ANALYSIS OF PROPOSED PREPROCESSING TECHNIQUE 

The block diagram of the basic preprocessing system is 
shown in Figure 1. The fundcimental design objective is to 
choose h (v) so that the functional 

00 

F = I u)(v)g^ [v)dv (1) 

— 00 

is minimized, where the "bar" over a variable indicates that 
the variable is a two-dimensional spatial vector and where 

oo 

g(v) = I 

— 00 

where denotes a convolution, subject to the constraints 

00 

Ki = 1 g^(v)dv (3) 

■ — 00 , • 

00 

K 2 = s(v)l^^(v)dv (4) 

'-.00 . - . 

K 3 = E {nT^(v)}_ (5) 

AS Stated previously, w(v) is an arbitrary penalty func- 

tion designed to influence the solution for h^ (v) so that 
g(v) is duration limited. The more rapidly w (v) increases 
with increasing v, the more rapidly g(v) will decrease with 
increasing v. Thus by choosing w(v) , g (v) can be made arbitrarily 
close to the desired impulse function. Similarly s (v) is an 
arbitrary penalty function designed to duration limit h^ (v) . 

Since any preprocessing will be performed digitally,' it is 
desirable to duration limit h^(v) so that truncation errors 
are minimized. 



- 9 - 


As a criterion for the degree of resolution improvement 
provided by a particular h^(v), the effective radius of the 
scanner aperture is defined as 


R ^ 
r' 


[ |v|^g^(v) 


dv 


nl 

7 


I 


dv 


( 6 ) 


Without preprocessing, the effective radius of the sceuiner 
aperture will be similarly defined as 




I |v|^h^(v)dv 




dv 


1 

7 


(7) 


Lagrange multipliers and the methods of functional analysii^'®^ 
will be used to solve Eq. 1 subject to the constraints of Eg. 3 
to 5. Eg. 1, 3, 4, and 5 may be combined into an augmented func- 
tional, I, which must be minimized with respect to h^(v) , 

00 00 00 

1 = 1 (i)(v)g^(v)dv + g^(v)dv + s(v)h^(v)dv 

,*00 *00 »oo 

+ A3E|n^(;)j 

00 00 00 ’ 

- I “(v) I h^(i)hjj(v - z)dz| hj.(u)hj^(v - G)dudv 

^ioo *00 — 00 • • 

00 «> 

+^lj I hj.(u)hjj(v - u)dudv 


—00 

00 


+A 2 I j s(u)6(u - v) h^ (v) hj. (u) dudv 
— 00 

00 . 00 

+Aje| I n(v - z)hj.(z)dz| n(v - u)hj.(u)du| 


( 8 ) 
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Eg. 8 'may be written in quadratic functional form as, 


1 = 11 aj^(u,z)hj.(z)hj.(u)dzdu 
*00 

00 


+Ai| 1 a2(u,z)h^(z)h^(u)dzdu 
*00 


00 

+A 2 I 1 a2(u,z)h^(z)h^(u)dzdu 
*00 


00 

*00 

(9) 

where aj^(u,z) , a 2 (u»z), a 2 (u,z] and a^(u,z) are 

linear opera- 

tors defined as°^ 

00 . 


aj^(5,5) = 1 a)(v)hjj(v - 5)hj^(v - 5)dv 
— 00 

(10) 

00 

a 2 (a,z) = 1 hjj(v - z)hjj(v - u)dv 
*00 '’ 

(11) 

1 - ■ ■ 

a 2 (u,z) a s(u)6(u - z) 

(12) 

a 4 (u,z) = E|n(v - z)n(v - u)| 

(13) 

= (z - u} , for n ( * ) a 



stationary ergodic random 
process. 

By taking the gradient of Eg. 9 with respect to 
where the adjoint linear operators of Eg. 10 to 13 are 

(14) 

(15) 


a2'(u,z) = 32 (u,z) 
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a2'(u,z) = a2(u,z) (16) 

and 

a^' (u,i) = a 4 (u,i) (17) 

and setting the gradient equal to zero, a homogeneous Fredholm 
integral equation of the second kind sometimes referred to as 
a Fredholm integral equation of the third kind, is obtained as 

00 

I [ai(v»z) + A^a2(v,S) + A3a4(v,z)jh^(z)d" 

■ «.00 • • 

+ A2s(v)hj.(v) = 0. (18)^^ 

The solution of this equation coupled with the constraint 
equations, Eq. 3, 4, and 5, would give the required point- 
spread function of the preprocessing filter, h^iv) . However, 
because of the numerical difficulties which may arise in the 
general solution of this type of equation and in order to more 
conveniently use the results of the multispectral scemner system 
analysis^ the solution for h^(v) will be formulated in the two- 
dimensional spatial frequency domain. 

(p) 

Eq. 9 may be rewritten using inner product notation as 

I “ + Aj^(a2h^,hj.) + A2(a3h^,hj.) 

+ A3(a4hr,hr) . (19) 

Z1 See Appendix A for derivation. 
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Eq. 19 may also be written as a quadratic functional in the 
spatial frequency domain as 

I = + ^2(®3”r'“r) 

+ ^3(B4H^,hJ (20) 

where H^(*) = 

^1*1 denotes the Fourier transfom 

and , B 2 / and B^ are the spatial frequency linear 

operators which are the Fourier transforms of the spatial 
domain linear operators a^^, a 2 » a.^, and a^. Thus 

f*f - . j2iTfu .2trva 

Bj^(f,v) =1 jaj^(u, 2 )e^ e^ dudz (21a) 

— oo • 

which after substituting Eq, 10 into Eq. 21a may be simplified 
to 

<21b) 

where ( • ) = ^ ^ } 

W(-) = 9'|u)(o| 

and Hj^*(») is the complex conjugate of Hj^(»). 

The adjoint of Bj^(f,v), defined as Bj^'(f,v), may be written as 

B^' (f ,v) = Bj^*(v,f) (22a) 
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which from Eq. 21b becomes 

= \(^}\*(f}w*(v - f ) . ( 22 b) 

Similarly, 



B 2 (f,v) = Hj^*(f)Hj^(v} 6 (f - 0 ) 

(23) 


B 2 *(f,v] = Hj^RH^*(f) 6 (v - f) 

(24) 


(f ,v) = S(f - v) 

(25) 

where 

S(«) = 3'^s(*) j , 



83 * (f ,v) = S*(v - f) 

(26) 


B^(f,v} = \*(v) <5(f “ v) 

(27) 

where 

♦ (•) = ?{r (•)} 



nn ■[ nn ' J 


and 


(28) 


The gradient of the quadratic functional of Eq. 20 becomes 


71 = (Bj^ + + A^(B 2 + 

+ A^fBj + + '‘3(24 + 84')^ 

which, upon expeinding the linear operator notation of Eq. 29a 
becomes 


^ 2 See Appendix B for a complete derivation of these 
spatial frequency linear operators. 
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’I - I [ + A^|b2(£,5) 

••00 

+ + A2|B3(f,C) + B3'(£,3)J 

+ A 3 /b^(£,v) + B^'(I,v)| ] H^(v)dv - 0. (29b) 

Substituting Eg. 21b, 22b, 23, 24, 25, 26, 27, and 28 into 
Eg. 29b, 

00 
— 00 

+ 1^2 {®(^ ■ ") + s*(v - I)|j H^(v)dv 

+ [jA^lHbCf)!^ + 2A3»nn(J)] Hr(J) - 0 . (30) 

Eg. 30 represents the general expression for the gradient 
of Eg. 20 with respect to which, when combined with the 

constraint eguations, completely specifies the spatial freguency 
spectrum of the preprocessing filter. The constraint eguations. 
Eg. 3, 4, and 5, may be rewritten in the spatial frequency domain 
as 

Ki = (31a) 

which after substituting Eg. 23 into 31a becomes 

00 

•'l " I 
-00 

00 

- J |Hb(?)l^|Hr(i)|^df, 


- v) dvdf 

(31b) 
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“ (BaHr'Hr) - (32a) 

which after substituting Eq. 25 into 32a becomes 

00 

^2 ” I " v)Hj.(v)H^*(f)dvdf (32b) 

••oo 

and 

K 3 - (33a) 

which after substituting Eq. 27 into 33a becomes 

00 • ■ • 

‘'3 “ n*nnM«(f - 5)H^(^)H^*(f)a5df , , 

— 00 

- ■ ■ 00 

“ [ • ; <”b>; 

— 00 ' 

Before Eq. 30 can be reduced to a form more suitable for 
the evaluation of the spatial frequency spectnam of the pre- 
processing filter, the penalty functions w(v) emd s(v) in 
Eq. 1 and 4 respectively, must be further examined. Since 
w(v) is designed to influence the solution of h^ (v) so that 
the composite imaging system point-spread function, g(v], is 
duration limited, ai possible choice for u)(v) would be 

(u(v) = 1 for £ V £ V 2 
= «>, otherwise. 

However, such a choice for w(v} would lead to analytical diffi- 
culties in Eq. 30, since the Fourier transform of aj(v) does not 
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exist. Thus the expression for w(v) must be chosen in such a 

manner that it allows enough flexibility to arbitrarily control 

the duration as well as the rate of decay of g(v] and, in addi- 

(28) 

tion, to have a Fourier transform' , Also, because of the form 
of Eg. 1, w(v} must be a positive valued function and also be 
convex to insure the existence of a global minimum to Eg. 1. 

One function for o)(v) which satisfies all the previous 

reguirements , written in terms of one variable, is 


for 



T2k 



0 < c < 1 


(34) 


and k a positive integer, as shown in Figure 2. For convenience, 
s(v) will also be described by the same type of function. 

The following analysis is based upon a rectangular coordinate 
system. Eg. 30 may be rewritten in terms of the x- and y- com- 
ponents of f and V, 

” = I j[ Hb*(*x'*y)"b(vV { - V’ ^ - V 

«-00 

+ W* (v^ - f„, - f ) > + A, { s(f - V , f - V ) 

X X y y' j 2 ( x x y y ■' 

+ S'K - ''y - *y)}] ”rK'V^''x'’'’y 

+ [ IVf'x'V * =0. 05) 
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For convenience in handling the analysis with respect to a 
rectangular coordinate system, both w (v) and s (v) will be 
defined as the product of their x- and y- components, from 
Eq. 34, 


‘^(v) = w„{x)w^(y) 


(36a) 


{ 


2x - - x^ 2 ,- -wx 


X ~ - X , 

w2 wl 


2k 


+ c 


( 


2y - y„i - y 


2k 


wl -^w2 1 'T 
— y + c. 


^w2 ^wl 


where 


0 < < 1 


0 < c < 1 

y 


(36b) 

(36c) 


and for k and k positive integers. 
Similarly, 


s(v) = s^(x)Sy(y) 
« 

{ 


(37a) 


• *sl ■ 


^s2 " ^sl 




+ d 


{ 


2k 


- ysi - ys2i =y 


ys2 - ^31 


(37b) 


where 


0 < d^ < 1 


0 < dy < 1 


(37c) 


and for kg^ and k^^ positive integers. 


Choosing k 


wx 


becomes 


1 = k „ , the Fourier transform of Eq. 36b 
wy ^ 



19 


W(f ,f ) = W [f )W [f ) 
'■ x' y' x'- X'' y '^ y' 


(38a) 


where 


«x(*x) = 


(x « - X , ) 
w2 wl 


1 

j 

TT 


2{x , + X 

6"(fJ - j — — 5'(fJ 


TT 


{Ki * ■" Kz - ’'wl^ ^x} ^(*x) 


(38b) 


and 


"y(^) = 


<y «2 - ywi> 


1 

tt' 


{..(, ) .j wl WZ_ J,(,j 


IT 


* ywz^^ * (ywz - y«i)^ =y} 


(38c) 


A similar expression results from taking the Fourier transform 
of Eq. 37b, 

Sff ,f 1 = S ff )S [f ) (39a) 

^ x' y' x'‘ X-' y'- y' 

where 
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and 


~ Z - y 


1 

■n‘ 


^ t^Sl ^82) 


«"(*y) - j ^ «-(fy) 


* {(^si ^ ^sa)^ ^ (^82 - ^Sl 




(39c) 


If it is assumed, as would usually be the case, that the penalty 
functions, a)(x,y) cuid s(x,y), are centered about the origin, 
then 


*wl ”^w2 ~ 
^wl ^w2 ~ 
’'si ”“’'s2 “ 

Ygi »-y32 “ 




Ys 



(39d) 


Substituting Eq. 39d into 38b, 

1 


■= - 


■ ”2 ■ 2 
4TT X 

W 


W (f ) = - 

y^ y; 


7*2 S’ 

4ir Y 


=x(*x) ■ - 
Sy (fy) - - 


7~T~7 
4ir X 

s 


“5™ 2’ 
y a 


38c, 39b, and 39c 

«"(fj + s«(fj 
«"(fj + 

«"(fy) + ay«(*y) 


f 


(39e) 

(39f) 

(39g) 

(39h) 
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Substituting Eq. 39e, 39f/ 39g, 39h, 39a, and 38a into Eq. 35, 


VI = 


w 


:a' 




■m T H. (f ,f Vh (f ,f ) 

2^i, 2 b*- x' y-' r'‘X' y-' 


- 4iT^y >r H. (f ,f )h ff ,f ) 

y aT 2 b '■ X ' y' r '■ x ' 


' w y 


3f 


y r X ' y’ 


- 4v^x *■ H. (f ,f )h ff ,f 1 

w X 2^ b'- X' y' r'^ x' y' 


3f 




-.4 2 ' 2 




>iw2 2, 

- 4tt X d. 

S X 


-*■ H (f ,f ) + 8 tt^x ^y ^ d d H ff ,f 1 
T r'“ x' y' s x y x' y' 


3f 


* [2^1 |Hb(f^,fy)l' ^ 2A34„„(f^,fy) 

which, when expanded in terms of H ff„,f }, becomes 

r X y 


Hr(*x-^y) ■ ® 
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r 1 ,f ] 


8 f 3 f 
X y 


+ 2A(£x-£y) 


+ 2A(fx'*y) 


3 H. (f ,f ) 
b X ' y-' 

3 ^H„(f ,f ) 
r». X- y' 

3 f 

3 f ^ 3 f 

y 

X y 

’“b(*x'*y) 

S^HrC^x-fy) 

3 f 

3 f 3 f ^ 

X 

X y 

3'Hb(£x,£^) 

3^Hr(£x,£y) 


3 f 3 f 
X y 


3 f 3 f 
X y 


3 


3 f. 


3 ^Hy.ff ,f ) 
r'^ x' v* 


3 f 


' *y) — (^X '*y) «b (*x ' *y) + ^ (^v -f„) 


3 f 


X ' y 2 


3 ^H (f ,f ) 
r X ' V'' 


3 f 


3^H, (f ,f ) — vi- /- 

2A(f^,£ J -y ■ ^ + 2E(f^,f J ^ y 

* y 3f 3f 
X y 


3 f 


3 H^ff ,f ] 
r '' x' 

3 f 


3 «K(f 


3 «H(f .fj 


2 A(f . 20 (f ,f 


^ y 3f 3f 

X y 


X' y' 


3 f 


3 H^(f »f ) 

r'^ x' V' 


3 f 


c(£x-£v) + 

^ y 3f„ 3f 


3 ^H.(f ,fj ,f ) 


X y 

> 2 , 


3 f 


X 


3 ‘^H. (f ,f ) ; 

+ E(f ,f } ■ . Ai . J £ ■ + H(f ,f ) 

x' y-' o '• x' y-' 


y 2 


“ 0 ( 40 b) 



23 


where 


H. ,f ) 

A(f f j . y 

w -^w 


(40c) 


B 


A^H. ,f ) 
rj: ^ V _ 2 b ^ x' y-' 

ifx'V " '“4^2 


8it^x y 
s -^s 


(40d) 


C(f ,f } = 2A. |h. (f ,f )I^ + 2A-4 (f ) (40e) 

'• X ' V'' 1 ' b X ' y' ‘ 3 nn X V'' 


)c 


2x‘x. 


w 


(40f) 


= - 


t:^ 


2" y, 


w 


(40g) 


A,d 

F(fx-fy) = - 

r 2x x 

s 


(40h) 


A»d 

= (*X'V = - 

^ 2ir 


(40i) 


Hff ,f 1 = c c |h, (f ,f ) 1^ + A„d d . (40 j) 

'• X Y'' X y' X y-' ' 2 x y 

Thus Eq. 40 in conjunction with the constraint equations, 

Eq. 31b, 32b, and 33b, specify the general form of the required 
preprocessing filter spatial frequency transform. 
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Separable Aperture 

A large class of physically realizable apertures 
may be modelled as separable apertures, where it is 
assumed that 


«b(*x’V = «bx(*x)»by(V 

(41) 

♦nn(*x'V = *nnx('x)*nny(«y) 

(42) 


(43) 


With the assumptions of Eg. 41 to 43, the solution of 

Eq. 40 can be considerably simplified by use of the 

method of separation of variables. Instead of substi- 

stuting Eq. 41, 42, and 43 into Eq. 40 and separating 

Eq. 40 into two differential equations, one a function 

of f cuid the other a function of f , a somewhat more 
* . Y 

fundamental approach will be used. 

Taking the inverse two-dimensional Fourier transform 

of Eq. 41, 

hb(x,y) = l^bx^^^^by 

Similarly, Eq. 42 and 43 become respectively, 

R „(x,y) = (y) (45) 

nn nnx nny 

and 

hj.(x,y) * hj.^(x)hj.y (y) . (46) 
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Writing Eq. 1 

in teras of the two spatial dimensions, 

oo . 

F » 1 |w(x,y)g^ (x,y)dxdy 
—00 

(47) 

where from Eq. 



g(x,y) = hj.(x,y)**hj^(x,y), 

(48) 

Substituting Eq. 44 and 46 into 48 and using the 
properties of two-dimensional convolution 

g(x,y) = 

“rx <’‘>*'ry <y>] ** [»bx <*> V 


a 

."rx <*> *\x <’'>] [*>ry <y> * V 'S'* ] 



(49a) 

where 

^x '*>**’«<*> *V'*> 

(49b) 


9y(y) = h^(y)*h^y(y). 

(49c) 


Substituting Eq. 36a and 49a intb' Eq. 47, 

OP 00 

— 00 —00 

» F F (50a) 

X y , ,,. ;. ;. .. , 
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where 


w 

= 1 “x 


(x)g (x)dx 

A 


% 

Fy = I Wy(y)gy(y)dy 
~00 


Substituting Eq. 49a into Eq. 3, 


= I gjj^ (x)dx I gy^(y)dy 
.00 .00 


= K, K, 
lx ly 


where 


lx 


w 

= 1 9.' 


(x)dx 


and 


ly 


I gy^ (y)dy . 


Substituting Eq. 37a and 46 into Eq. 4, 


00 oo 

.oo .oo 


= ?'2x^2y 


where 


w 

1 


^2x “ 


and 


K 


2y 


w 

I =y(y»v. 


(y)dy 


(50b) 

(50c) 

(51a) 

(51b) 

(51c) 


(52a) 

(52b) 


(52c) 
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Substituting Eg. 42 and 43 into Eg. 33b,, : 


00 . oo 

‘'3 = 1 »nnx(*x) l «rx (*x) -l ^«x f ♦nnyCV 


where K, = { 4 (f )|h ff )|^ df 
3x J nnx'- x-' ' rx*- x' ‘ x 

00 

00 

= I 4 „ (f )df 
J nnTx'- x' X 


(53a) 

(53b) 

(53c) 


where *nnTx^^x^ power spectral density of the 

x-component noise in the processed image, and since 
the noise is assumed to be a sample function from a 
stationary ergodic random process, 

•'3X = '’{"tx"<=‘>}. 

Similarly, 

00 

■'3y - I "y 

.00 

- E{"iy^(y)}_ 

Thus , the problem of determining the optimum prepro- 
cessing filter point-spread function, h^ (x,y) , reduces 
to finding the h (x) that will minimize Eg. 50b 


(54a) 


(54b) 
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subject to the constraints of Eq. 51b, 52b, and 53d, 
and to finding the h^y(y) that will minimize Eq. 50c 
subject to the constraints of Eq. 51c, 52c, and 54b. 
The original two-dimensional preprocessing problem 
reduces to two one-dimensional processes which have 
similar equations. From the proceeding vector nota- 
tional analysis used in Eq. 1 to 34 and Eq. 36 to 39, 
the system of equations necessary to solve for h (x) 

i A 

and bj.y(y) may be formulated. 

To solve for h_ (x) , the augmented quadratic 
functional of the form of Eq. 8 determined by Eq. 50b, 
51b, 52b, and 53d becomes 

OP 00 

{ * hxf 

*00 ••oo 

00 

«00 

which may be written as a quadratic functional in 
the spatial frequency domain from Eq. 20 as, 

^x - (^ix^rx'^rx) ^ ''lx(®2x"rx'«rx) 

* ''2x(B3x»rx'«rx) + ^3x(»4x«rx'“rx) 
where from Eq. 21b, 

hx(^x'''x) - Ax'f^xKxKJWxt^x - '’x) 



(57) 
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from 

Eq.. 

22b, 

. : • 'T • 



B' 

lx 


= »bx(''J«bx*(*x)«x*(''x - *x) 

(58) 

from 

Eq. 

23, 




®2x 

(*x'''x) 

" “bx*(*x)»bx('’x)«(«x - '-x) 

(59) 

from 

Eq. 

24, 




®ix 


- «bx(''x)«bx*(£x)«(''x - *x) 

(60) 

from 

Eq. 

25, 





«3x 

= ®x(‘x - '’x). 

(61) 

from 

Eq. 

26 , 





Bl 

3x 

(V-’x) = ®x*('’x - *x) ’ 

(62) 

from 

Eq. 

21 , 





^x(‘x' 

x-* nnx X' '•X x' 

(63) 

and from 

Eq. 28 , 





B4x(^x'' 

^ ) = * ff ]6fv - f 1 . 

x-* nnx'^ X' ^ X x-' 

(64) 
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The gradient of Eq. 56 may be written in the 
form of Eq. 29a 

- Kx * «lx')»rx ^lx(®2x ■^>2x')»rx 

'^2x(®3x ®3x'^“rx ^ ''3x(«4x «4x’^“rx ' <«5> 

which may be expeuided to the form of Eq. 30, 

.00 

’^x = 1 [«bx* (*x) "bx (%) {«x (*x - ''x^ + ”x* (“x - *x} 

. , 00 . 

+ '^2x {®x(*x - ''xi + Sx'^x - ^x)}]»rx('’x)^% 

* [2''lxl“bx(*x)l' * = 0. (66) 


Substituting Eq. 39e and 39g into Eq. 66, the following 
differential equation arises 



* 0. (67p 


/,3 For derivation see Appendix C. 



From Eq. 31 and 51, 


w- 

■'lx ” l.l■^bx(yl'l■•rx(^x)l' “x <««> 


cuid from Eq. 32, 39g, and 52 


K 


2x 


4tt X 


CO 

5^1 [“;rx(*x)'‘rrx(*x) 


s —CO 


+ H”. (f )h . [f )ldf 
rix^ x-' rxx'^ x-'J x 

00 

+ <*x I l"rx(*x)l^ «x, 

■■ — » 

where + JHrixt^x) . 

and where H (f ) is the real part and H . (f ) is 
iTJTx X irxx X 

the imaginary part of 

Restating Eq. 53b, 

00 ■ 

■'3X ' j ♦nnx(*x)l«rx(*x^ l^ ^*x • 

*00 

Thus the simultaneous solution of the differential 
equation Eq. 67, and the constraint equations, 

Eq. 68, 69, and 53b, specify the form of the x-ccanpon- 
ent of the spatial frequency transform, or equivalently 


Y4 For derivation see Appendix D. 
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the point-spread function, of the preprocessing 
filter. 

In a similar manner it is possible to solve 
for h^y(y) by forming an augmented quadratic func- 
tional of the form of Eq. 8 determined by Eq. 50c, 51c 
52c, and 54b, 

CO CO 

= j ‘*»y(y)gy^ (y)dy + ^ly j ^y^fyJ^^y 

— 00 «,oo 

' oo'-, 

Vi Sy<y>*'ry^'y>‘3y + "sy^Ky^'y’}- '’®> 


By following an analogous procedure to that used for 
determining h „(x) in Eq. 55-69, the equations which 
specify h^y (y) may be formulated. Only the results 
will be stated since the derivation of the equations 
for h^y(y) is identical in form to that given for 
h^^(x) with the appropriate change in variables from 
X- to y- dependency . 

The differential equation specifying the form 
of H^y(fy) becomes. 


+ y 


H" (f J + 5 H (f ) 

[«by* (fyKy (fy) (°y^''ly) I “by (fy) > ^*''3y*nny(^y> «y}J 


s 


•H 


l“by(*y)l^ + ''2yV 


CfJ 


ry^ y 


= 0 




(71) 
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while the constraint equations become, 


w 

•'ly - ] 


( 72 ) 


K 


2y 


4tt y 


00 

^ i bv (^y) ^rry ^^y^ 


OD 

h”. ff )h . (f M df + d I |h (f )|^ df 

rxy'> y' rxy'- y-'J y y J ry'- y-' ' i 


and K- = 
3y 


$ (f ) |h (f ) r df . 
nny' y' ' ry'^ y-' ' y 


(54a) 


Thus the simultaneous solution of Eq. 71, 72, 73, and 
54a will specify the y-component of the preprocessing 
filter. 

2 , Radially Symmetric Aperture 

Probably the most common type of aperture, because 
of the physical ease in construction, is the radially 
symmetric aperture. For this case it is assumed that 


H, If ,f 1 = 

b'' x' y' 


(74) 

(f ,f ) = 
nn'' x' 

♦nnr(*r) 

(75) 


H (f ) 
rr'- r' 

(76) 

f ^ + f ^ = 

X y 


(77) 


(73) 


where 
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The solution to Eq. 40 is analogous to the 
solution for where all variables dependent 

upon X are replaced by corresponding variables 
dependent upon r. From Eq. 36a, 

w(v) = o)^(r) 

, r - ’^w2 

^w2 ^wl 


for a positive integer auid 




2 

r = 

2^2 
X + y . 


(78b) 

From Eq. 

37a, 






s (vj * s^ (r) 





3 

2r - 

^sl - ^s2 ' 

2k 

sr 

(79) 



L ^^s2 - ^sl J 



for 

k^^ a positive 

integer. 



Choosing 


k = 
wr 

1 “ kgj. and 




^wl = 

-r = r 

w2 w 


(80) 


^sl = -^s2 = ^s' 


-I 2k 


wr 


(78a) 
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the Fourier transforms of. Eq.. 78a and 79 become 


W (f ) « — =• 6'’ff ) 

w 


( 81 ) 




(82) 


4ir^r 


To solve for (r ) , the augmented quadratic 
functional of the form of Eq. 8 which must be mini- 
mized with respect to h (r) becomes, 

.jT-i • • •’ . 

OO 00 

Ir “ j ^Ir | 

00 

+ j =^(r)h^/ (r)ar 

— 00 

* '^3r ^ . <«) 

By following an analogous procedure to that used 
for determining h (x) in Eq. 55-69, the equations 

Jl X 

which specify h^^(r) may be formulated. Again, only 
the results will be stated since the derivation of 
the equations for h^^(r) is identical in form to that 
given for h (x) with the appropriate change in 
variables. 

The differential equation' specifying the form 
of H^^(f^) becomes. 
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V . NUMERICAL SOLUTION TECHNIQUE FOR PREPROCESSING; FILTER 


Since the form of the differential equations specifying 
the shape of ' \y(v) r Eq. 71, and h^^(r), Eq. 83, 

is the same, a single method of solution is applicable. In this 
discussion specific reference will be made to the solution of 
Eq. 67, 68, 69, and 53b for H (f ) . 

Eq. 67 may be written in the form. 


H" (f ) + A(f )h' (f ) + B(f )h ff } = 0 


X-' rx'- X 


X-' rx*- x-* 


(87) 


where 


s 


|H, (f )l^x ^ + A, X ^ 
' bx '■ x-' ' s 2x w 


( 88 ) 


B(fx) 


and 

, 2 [Hbx*(^x)“bx(^x)-^'''’'w’f (V*lx) l”bxt^x) I ’'•'^ 3 x'*'nnx(*x)''^x}] 


«bx(*x)l“^ 


+ A- X ' 

2x w 


It should be noted that H (f ) , A(f ) and B(f ) are complex. 

ITX X X X 

Defining 


"rxtfx) = «rrx(*x) , 

; (69b) 

*(*x) - *r(*x) jM^x) . 

(90) 

and 



'(91) 

as the sum of real and imaginary components , then 
Eq. 69b, 90, arid 91 into Eq. 87, 

by substituting 

[»;rx(*x) ^ + [a^(£^) + 

'■ [®r(y + ^®i(^x)][«rrx(^x^ 2H^ix(y_ 

=0 (92) 


(89) 
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or by collecting the real and imaginary components, 


H" ff ) + A (f )h' (f^l - A. (f )h*. (f ) + B (f )h (f ) 
rrx'- x-' r*- x-' rrx'“ x' i '■ x-' rix'^ x-' r*- x' rrx'- x-' 

- B. (f )h . (f ) + j|H"- ff ) + A. ff )h' ff ) 

+ A ff )h’ . ff ) + B. ff )h ff ) + B ff )h . ff ) 
r'- X-' rix'- x' 1 x' rrx'- x-' r'- x-' rix'- x' 

(93) 

Eq. 93 may be separated into two differential equations formed 

(17) 

by the real and imaginary components of Eq. 93 , 


= 0. 


H" If ) + A ff )h' ff ) - A. (f )h' . ff ) + B ff )h ff ) 
rrx'- x-' r*- x-' rrx'- x-* x'- x-' rxx> r'- x-* rrx*- x-' 




(94) 


and 


H". ff } + A. ff )h' ff ) + A ff )h' . ff 1 + B. ff )h ff ) 
rix'- X-' x'- X-' rrx'- x' r'- x' rxx'- x-' x'- x' rrx'- x' 


B ff )h . ff } = 0, 
r'- X-' rxx'- x' 


(95) 


Thus, the original complex second order differential equation, 

Eq. 87, has been reduced to a system of second order differential 
equations, Eq. 94, 95. 

To make use of the many subprograms available for handling 
systems of first order differential equations, Eq. 94 and 95 may 
be reduced to a system of first-order differential equations 
by introducing the variables, 


«l(^x) V«rrx(«x) 
«2(fx) = Hrrx(*x) 

= »rix(*x) 

«4(y = »kx(*x) 
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By substituting Eq. 96 into 94 and 95, the following system of 
first-order differential equations is formed. • 





■>2 <*x' = - ‘>1'<*X>. - \'*x> W2 «x> + »i «x>^ »3 '*x> 


+ A. (f^) (f^) 

«3<*x> = «4'V 

HJ(fx) = - %<*x> 

- ^r'*x> «4<*x> 


A. (f ) H, (f ) 

I X 2 X 


B (f ) (,f ) 

r X 3 X 


(97) 


After specifying initial conditions for Hj^(O), H 2 (O) , ^^(O) 
and H. (O) , H (f ) may be obtained as 

4 rx X 


Hrx'*x> = »l<fx> 3 “3<*x> 


(98) 


In order to solve for H (f) , the system of differential 
equations, Eq. 97, plus the constraint equations, Eqw 67, 69, 
and 53b, must be solved simultaneously. The constraint equations 
may be considered to be a system of non-linear equations where 
the unknown parameters are ^ 2 x' ^3x* For a given value 

of these parameters, Eq.97 may be used to determine H (f ) 
and the constraint equations checked to determine if they are 
satisfied. If the constraint equations are not satisfied, 
appropriate pertubations in ^ 2 x' ^ 3 x made and a 

new value of H(f) computed. This procedure would be repeated 

ITX X 

until the constraint equations are satisfied. A program for 
solving a system of non-linear equations has been developed and 
could be used for determining ^ 3 x* 

Two possible problems which might prevent obtaining a solu- 
tion for H (f ) using the procedure described above must be 
considered. One problem would be the possibility of obtaining 
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a solution which would represent a local minimum to Eq. 55 but 
not necessarily the best solution which would represent a 
global minimum. This problem should not arise because of the 
choice of the functionals in Eq. 50b, 51b, 52b, and 53d. Since 
all of these functionals are convex, a global minimum is assured 
The second problem that could arise is related to the 
solution of the system of differential equations, Eq. 97. Is 


(16) 


it possible for a given set of ^2x' 


and A, that either no 
3x 


solution or several solutions to Eq. 97 exists? 
Section 7 of 


A theorem in 


states that as long as B_(f ), B. (f ) , 




and continuous functions, then for a given set of 

initial conditions one cind only one set of solutions exist. 



VI. CONCLUSIONS. .. , ,, 

The preceeding analysis was based upon k = k = k = k = 1 , 
or for ^sr~ If a convenience for formulating the equations 

specifying the required preprocessing filter. The larger the 
integer value of k^ or k^, the greater the reduction in the effec- 
tive scanner radius of the composite imagery system. However, 
the size of the system of differential equations increases and 

results in a system of 4k linear differential equations , where 

max 

k is the largest integer value of k k , k , k , or k , 
max ^ ^ wx , wy ' sx ' sy ' wr ' 

k . For example, ifk =3,k =2,k =k =1, then a system of 

sr wx wy sx 

differential equations similar to equation 97 will result composed 

of 12 simultaneous linear differential equations. However, 

one of the major advantages of this proposed technique for image, 

preprocessing is that it is sufficiently general to allow for 

any integer value of k and k and thus allows the effective 

w s 

scanner radius of the composite imaging system to be arbitrarily 
reduced, subject primarily to the noise constraint. 

Although Equation 40 in conjunction with Equation 31b, 32b, 
and 33b, specify the general form of the spatial frequency trans- 
form of cuiy preprocessing filter, it was shown that in the case 
of a separable aperture or a radially symmetric aperture the solu- 
tion can be considerable simplified. These two classes represent 
the most common types of apertures used for data collection. 

There exist many functions which are separable in the sense 
that the function can be expressed as the product of its x- and 
y- components. By appropriately choosing such a function, many 
types of asymmetric apertures may be approximated as symmetric, 
separable apertures. For example, by properly selecting the para- 
meters of a two-dimensional Gaussian function, an elliptical aper- 
ture of uniform density could be approximated. By using a two- 
dimensional Gaussian function, it is also possible to approximate 
a radially symmetric aperture as a separable aperture. The 



42 


principal advantage for choosing a separeible aperture is that 
two-dimensional convolutions with such an aperture is equivalent 
to two one-dimensional convolutions along each orthogonal axis. 
Thus preprocessing time for a given data set can be significantly 
reduced by using a separeible aperture. 

The preprocessing filter theory presented in this paper, 
and in particular the filters described in Equations 67, 71, and 
83, will be applied to the multispectral scanner data at LARS to 
determine the best set of parameters for reducing the effective 
scanner aperture and the effect of such a reduction on classifi- 
cation accuracy. 



43 


References 


1. Andrews, H. C. , A. G. Tescher, and R. P. Kruger, "linage 

Processing by Digital Computer" , I.E.E.E. Spectrum , 
July 1972, pp. 20-32. 

2. Arguello, R^ , H. Sellner, and J.. A. Stuller, "Transfer 

Function Compensation of Sampled Imagery", I.E.E.E. 
Transactions on Computers , July 1972, pp. 812-818. 

3. Becherer, R. J. , "Optimum Shaded Apertures for Reducing 

Photographic Grain Noise", SPIE - Image Information 
Recovery Seminar Proceedings, Oct. 54-25 , 1968, pp. 
■ST-96. 

4. Bracewell, R. N. , "Correcting for Gaussian Aerial Snooth- 

ing" , Austrialian Journal of Physics , 8:54, 1955, 
pp. 54-66. 

5. Bracewell, R. N. , "A Method of Correcting the Broadening 

of X-ray Line Profiles", Austrialian Journal of 
Physics , 8:61, 1955, pp. 61-67. ~~ 

6. Bracewell, R. N., "Correcting for Running Means by Succes- 

sive Substitutions", Austrialian Journal of Physics , 
Vol. 8, No. 3, 1955, pp. 329-334 . "" 

7. Epstein, B. , Linear Functional Analysis, W. B. , Saunders, 

Co. , 1970. 

8. Franks, L. E. , Signal Theory , Prentice-Hall, 1969. 

9. Harris, J. L. Sr., "Image Evaluation and Restoration", 

Journal of the Optical Society of T^erica ", May 1966, 
pp. 569-574 . ’ 

10. Harris, J. L. Sr. , "image Restoration", Proceedings of the 

National Aero^ace Electronics Conference , 1963 , pp. 

WT-m . — — — ^ 

11. Horner, J. L. and M. J. Rycus, "Restoration of Linearly 

Smeared Transparencies", SPIE - Image, Inf oration 
Recovery Seminar Proceedings, Oct. 24-25, 196^8 , pp. 
TIT- 129. — 

12. Huang, T. S., W. F. Schreiber, et. al. , "Image Processing" 

Proceedings of the I.E.E.E., Nov. 1971, pp. 1586-1609 



44 - 


13. Kinzly, R. E., R. C. Haas, and P. G. Roetling, "Designing 

Filters for Image Processing to Recover Detail", SPIE - 
Image Infor matio n Recov ery Seminar Proceedings, Oct. 

24-25 , 1968 , pp. 97-104 . 

14. McGillem, C. D. and T. E. Riemer, "Moire^ Patterns and Tv/o- 

Dimensional Aliasing in Line Scanner Data Acquisition 
Systems", LARS Information Mote 111772 . 

15. Nathan, R. , R. Selzer, and F. Billingsley, "Image Processing", 

JPL Technical Report 30-792 , JanlO, 1970. 

16. Pierre, D. A., Optimization Theory with Applications, Wiley, 

1969. ■ ^ 

17. Pontryagin, L. S. , Ordinary Differential Equations , Add i son - 

Wesley, 1962. 

18. Pratt, W. K. , "Generalized Wiener Filtering Computation 

Techniques", I.E.E.E. Transactions on Computers, July 
1972, pp. 636-T41. 

19. Rice, R. B. , "Inverse Convolution Filters", Geophysics, Vol. 

XXVII, No. 1, Feb. 1962, pp. 4-18. 

20. Robbins, G. M. and T. S. Huang "Inverse Filtering for Linear 

Shift-Variant Imaging Systems", Proceedings of the I.E.E.E. , 
July 1972, pp. 862-872. 

21. Robinson, E. A., Statistical Communication and Detection , 

Charles Griffin and Co., 1^67. 

22. Slepian, D. , "Restoration of Photographs Blurred by Image 

Motion", Bell System Technical Journal, Dec. 1967, pp. 
2353-2362.. 

23. Sondhi, M. M. , "Image Restoration: The Removal of Spatially 

Invariant Degradations", Proceedings of the I.E.E.E. , 

July 1967, pp. 842-853. 

24. Smith, H. A., "Improvement of the Resolution of a Linear 

Scanning Device", SIAM Journal on Applied Mathematics , 

Jan. 1966, Vol. 14 , No. 1. ^ ^ 

25. Stuller, J. A., "An Algebraic Approach to Image Restoration 

Filter Design" , Computer Graphics and Image Processing , 

Dec. 1972, pp. 107-122. ^ 

26. Yosida, K. , Lectures on Differential arid Integral Equations 

Interscience Publishers, Inc. , 1960 . 



45 


27. Yu, F. T. S., "Image Restoration, Uncertainty, and Infor- 

mation", Applied Optics , Jan, 1969, pp. 53-58. 

28. Zemanian, A. H., Distribution Theory and Transform Analysis , 

McGraw-Hill, 1965. 


Bibliography 

1. Dettman, J. W. , Mathematical Methods in Physics and Engineering, 

McGraw-Hill, 1962. ^ ^ 

2. Papoulis, A., The Fourier Integral and Its Applications , 

McGraw-Hill , l$6i . 

3. Papoulis, A., Systems ^d Transforms with Applications in 

Optics , McGraw-Hill , 1968. 

4. Temes, G. A,, V. Barcilon, and F. C. Marshall III, The Optimi- 

zation of Bandlimited Systems", Proceedings of I.E.E.E. , 
Feb 1973, pp. 196-234. 



46 


APPENDIX A - DERIVATION OF EQUATION 18 

The gradient of Eq. 9 may be written as 
00 

VI = f (a^^ (u, z) + aj^ {u, z) ] hj. (z) dz 


—00 

00 


+ 

''i ; 

[a 2 (u, z) 

+ 

*2 

(u, z)] 

hr 

(z) 

dz 



— 0 
0 

0 

0 

f 








+ 

^2 

' [a- (u, z) 

1 - J , 

-t- 

^3 

(u, z)] 

hr 

(z) 

dz 



— 0 
0 

O' 

0 

f 








+ 

''3 . 

— 0 

[a^ (u, z) 

O'' 

-f 

*4 

(u, z) ] 

hr 

(z) 

dz . 

(Al) 


Substituting Eq.- 12, l4, 15, 16, and 17 with Eq. Al, 

00 

VI = j (u, z) + ^1-^2 ^3 ®4 ^r 

*“00 

oo 

+ A 2 I s(u) 6(u - z) (z) dz 

— 00 
00 
* ' 

= [a^^ (u, z) + (u, z) + (u, z) ] h^ (z) dz 

— CO , 

+ A 2 s(u) hj,(u) = 0 . (A2) 
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APPENDIX B - DERIVATION OF EQUATIONS 21b, 
23, 25, and 27 

Siabstituting Eq. 10 into 21a, 

00 

Bj^(f,v) “II w(v) (v - i) hj^ (v - u) dv gj27rVz 


Introducing a change of variable in Eq’ B1 'where ' 


V - u = a 
00 


» / 

B^(f,v) = 


w (v) h, (v - z) h, (a) a) ^j2irvz 

b b 


= (t) 


II w (v 


) hj^ (v - z) dv dz 


and introducing another chauige of variable in Eq. B3 where 

V - z = B 

00 

I viv) hjj(6) dJ dS 


= H^( 


v) w(v 


-) -j2TTV(f - V) 


dv 


Following a similar development 
00 ' , 

B2(f, v) = I a2(u, z) du dz 


which after substituting Eq. 11 into B6a becomes 


B2 (f 


hj^ (v - z) hj^(v - u) dv du dz 


(Bl) 


(B2) 


(B3) 


(B4) 


, (B5) 


(B6a) 


(B6b) 
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Introducing the change of varicible of Eq. B2 into B6b, 


B2(f,v) = 


(I \ - 


z) h^^(a) -«) e^^TTVz 


= «b* 


(f) II (y ~ z) gj2nvz 


cind introducing the change of variable of Eq. B4 into B7, 


B 2 (f ,v) = Hj^*(f) 


h^(6) '6) av dS 




-j2i,v(( -V) 


= Hj^*(f) Hj^(v) 6'(f -V). 


B3(f , V) = 


a3(G, i) ij2Trfu ^j2TTVz ^ 


Substituting Eq. 12 into B9 , 


B3(f, V) 


= I S (u) 6 (u - z) gj^TTVZ 


s 


(5) i^2wi(f -^1) 


= s(f -V) .: 


B 


(f, U) =. II 3^(5', i) ;i2*K: , j?.vi ^5 


Substituting Eq. 13 into Bll, 


00 
r f 


B^(f, v) = 


R (5 - S) 532r£u d5 di . 

nn 


(B7) 


(B8) 

(B9) 


(BIO) 

(BID 

(B12) 


~.00 



4 9. - 


Introducing the change of. variable 


into Eg. B12 , 



00 
/ f 


a^(f, V) 


^ <35 d5 

nn 


= <t * (v) 
nn 


5j2itu(f -V) 


= « *(v) 6(f - 5). 

nn 


(B13), 


(B14) 
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APPENDIX C - DERIVATION OF EQUATION 67 
Sxibstituting Eg. 39e and 39g into Eq, 66 , 


VI 


= 1 “bx'V 


— 6" (f - V ) + c 6 (f - V ) 

2 X X X X X 


W 


6 " (v . - f ) + c 5 (v - ' f ) 

,22 X X' X X X 


+ A_ ( ^ 5" (f„ - v„) + d..6(f.. - v„) 

V Arr^ V 


2x 


47T*X 

s 


XX XXX 


47T X 


6" (V - f ) + d 6(v 
2„ 2 X X X X 


- V } 
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H (v ) dv 
rx X X 


2A, H, (f ) ^ + 2A, $ (f ) 

lx bx X 3x nnx x 


H (f ) = 0, 
rx X ' 


(Cl) 


or 


\x*(^x^ d^ 

\7I = — — ~ [H, (f ) H (f ) ] 

2tt^x df ^ ^ ^ 

W X 


+ 2c H, (f ) H (f ) 
X ' bx X ' rx X 


^2x d^ 

H (f ) + 2d H (f ) 

2tt^x ^ df ^ ^ X rx X 

S X 


* |«bx'vr' 2''3 x ♦nnx>'x>' “rx'^x* = 


(C2) 
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or, 

VI 


H, * ( f ) 
bx X 

27 t 2 x 2 

w 


[H, " (f ) H (f ) + 2H, ' (f ) H ^Xf .) . 

bx X rx X bx x rx x 


A -) ■ ' : r-‘ 

+ H, (f ) H ( f , ).] - - ~( f ) 

bx X rx X 2n^-x * ’ ' 

■ s 


■X ■! 


+ [2(c + A, ) H, (f ) |2 + 2A, $ (f ) + 2d ] H (f ) = 0, (C3) 

^ X lx bx X I 3x nnx x x rx x ' 
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or in normalized form, 

2H. * (f ) H, '(f) X ^ 

W = H " (f ) + bx X _. s_ 

|H^ (f )]' X “ + X, X ” * 

' bx X ' s 2x w 
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1 

(C6) 
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APPENDIX D - DERIVATION OF EQUATION 69 
From Eq. 32, 39g, and 52, 
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From Eq. 69b, 
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Substituting Eq. D5 into D2 , 
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Since h (x) is assumed to be a real function, then 
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H” . (f„) is an odd function, 
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Thus, the second integral in Eq. D6 is zero, and 
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